What is claimed is: 



1. A recursive discrete Fourier transformation method 
wherein dat a values x(t) , x(t+l) , x(t+2) , x(t+3) , . . . , x(t+N-l) , 
5 x{t+N) sampled at times t, t+1, t+2, t+3, t+N-1, t+N (N 

is a positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 
since time t as a data stream, a frequency component, which 
is degree k ( k is 0 or a positive integer smaller than N ) obtained 
10 by carrying out complex Fourier transformation on the data 
stream, is obtained such that a real part Xr(k, t) thereof 
and an imaginary part Xi(k, t) thereof are complex Fourier 
coefficients, the method comprising: 



15 x(t+2), x(t+3), x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients Xr{k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 

20 a third step of storing the complex Fourier coefficients 

Xr(k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using data value x{t+N) supplied 
at time t+N, data value x(t) stored in the first memory means 

25 temporarily and the complex Fourier coefficients Xr(k, t), 
Xi(k,t) stored in the second memory means temporarily, 
obtaining complex Fourier coefficients Xr{k, t+1) and Xi(k, 
t + 1) to the data stream supplied since time t + 1 with respect 
to a positive constant value A for giving an amplitude value 

30 corresponding to a difference between the x(t+N) and the x(t) 
according to a following equation: 



a first step of storing the data stream x(t), x(t+l). 




jik 



+l) = Z.(^,/)cos 2 
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-43- 



2. A recursive discrete Fourier transformation method 
according to claim 1 wherein inverse discrete Fourier 
transformation data subjected to inverse discrete Fourier 

5 transformation using a positive constant a is supplied so as 
to obtain a complex Fourier transformation coefficient to the 
supplied inverse discrete Fourier transformation data, 

the discrete Fourier transformation being carried out 
using the positive constant A corresponding to the positive 
10 constant a. 

3. A recursive discrete Fourier transformation method 
according to claim 1 wherein k indicating a degree of the complex 
Fourier coefficient is set to a desired value and Fourier 

15 transformation to the set value k is repeated at every time 
having an equal interval. 

4 . A recursive discrete Fourier transformation method 
wherein with sampling frequency as fs , data values x(t ) , x(t+l), 

20 x(t + 2) , x(t4-3) , . . . , x(t+N-l) , x(t+N) sampled at times t , t + 1, 
t+2, t + 3, t+N-1, t+N (N is a positive integer which is 

1 or more) each having an equal interval are supplied and with 
such N data values supplied since time t as a data stream, 
a frequency interval given with a minimum frequency fl and 

25 a maximum frequency f2 to the data stream is regarded as a 
measuring frequency interval, a frequency interval obtained 
by dividing that measuring frequency interval by the N is 
assumed to be an analysis frequency interval and a result of 
the frequency analysis obtained by carrying out complex Fourier 

30 transformation at every analysis frequency interval is 
obtained as a frequency component which is k ( k is 0 or a positive 
integer smaller than N) times the frequency interval, such 
that a real part Xr(k, t) thereof and an imaginary part Xi(k, 
t) thereof are complex Fourier coefficients, the method 

35 comprising: 

a first step of storing the data stream x(t), x(t+l). 
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15 



x(t+2), x(t+3), x{t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients Xr{k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients 
Xr(k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using data value x(t+N) supplied 
at time t+N, data value x(t) stored in the first memory means 
temporarily and the complex Fourier coefficients Xr(k, t), 
Xi(k,t) stored in the second memory means temporarily, 
obtaining complex Fourier coefficients Xr(k:, t+1) and Xi(k, 
t+1) within the frequency interval with the minimum frequency 
f 1 and the maximum frequency f 2 to the data stream supplied 
since time t+1 with respect to a positive constant A for giving 
an amplitude corresponding to a difference between the x(t+N) 
and the x(t) according to following equations: 



20 



X^{k,t +l)=j^X^{k,t) + ~[x{t +N)~x{t)'}^xcos |2 



sin J 2 



{f2-fl)k 



N-1 



+ /1 



if2-fl)k 



N-1 



+ /1 



Z,. + 1) = Z, (A;,? )x cos J 2— 



if2-n)k 



N-1 



+ /1 



X,{k,t) + ^[x{t + N)-x{t)] 



xsm J2 — 



(/2-/i)A: 



+ /1 



5. A recursive discrete Fourier transformation method 
25 according to claim 4 wherein inverse discrete Fourier 
transformation data subjected to inverse discrete Fourier 
transformation using a positive constant a is supplied so as 
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to obtain a complex Fourier transformation coefficient to the 
supplied inverse discrete Fourier transformation data, 

the discrete Fourier transformation being carried out 
using the positive constant A corresponding to the positive 
5 constant a. 

6. A recursive discrete Fourier transformation method 
according to claim 4 wherein k indicating a degree of the complex 
Fourier coefficient is set to a desired value and Fourier 

10 transformation to the set value k is repeated at every time 
having an equal interval. 

7. A recursive discrete Fourier transformation method 
II wherein data values x(t) , x(t+l) , x(t+2) , x(t + 3) , . . . , x(t+N-l) , 
^ 15 x(t+N) sampled at times t, t+1, t+2, t+3, t+N-1, t+N (N 

is a positive integer which is 1 or more) each having an equal 
-;f interval are supplied and with such N data values supplied 

since time t as a data stream, a frequency component, which 
is degree k (k is 0 or a positive integer smaller thanN) obtained 
20 by carrying out complex Fourier transformation on the data 
stream, is obtained such that a real part Xr(k, t) thereof 
and an imaginary part Xi(k, t) thereof are complex Fourier 
coefficients, the method comprising: 

a first step of storing the data stream x(t), x(t+l), 
25 x(t+2), x(t+3), x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients Xr(k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 
30 a third step of storing the complex Fourier coefficients 

Xr(k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the complex Fourier 
coefficients Xr(k, t) - jXi(k, t) stored in the second memory 
35 means temporarily, obtaining complex Fourier coefficients 
Xr(k, t+1) - jXi(k, t+1) to the data stream supplied since 
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time t+1 based on a transfer function expressed in a following 
equation , 



/f(z) = ^(l-z-^) 

where A is a positive constant for providing [x(t+N) - x(t)] 
5 with an amplitude. 

8 . A recursive discrete Fourier transformation method 
according to claim 7 wherein the positive constant value A 
for providing with an amplitude corresponding to a difference 

^ 10 between the x(t+N) and the x(t) is capable of being set 

CI, selectively with an inverse number of square root of 1 , N or 

W 1/N. 

9, A recursive discrete Fourier transformation method 
g 15 according to claim 7 wherein k indicating the degree of the 
ij complex Fourier coefficient is set to a desired value and the 
'4 Fourier transformation to the set value k is repeated at every 
J time having an equal interval. 

20 10* A recursive discrete Fourier transformation method 

wherein with sampling frequency as fs, data values x{t ) , x{t + l) , 
x(t+2),x(t + 3), x(t+N-l) , x( t+N) sampled at times t , t + 1 , 

t + 2, t+3, t+N-1, t+N (N is a positive integer which is 

1 or more) each having an equal interval are supplied and with 

25 such N data values supplied since time t as a data stream, 
a frequency interval given with a minimum frequency fl and 
a maximum frequency f 2 to the data stream is regarded as a 
measuring frequency interval, a frequency interval obtained 
by dividing that measuring frequency interval by the N is 

30 assumed to be an analysis frequency interval and a result of 
the frequency analysis provided by carrying out complex Fourier 
transformation at every analysis frequency interval is 
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obtained as a frequency component which is k ( k is 0 or a positive 
integer smaller than N) times the frequency interval, such 
that a real part Xr(k, t) thereof and an imaginary part Xi(k, 
t) thereof are complex Fourier coefficients, the method 
comprising: 

a first step of storing the data stream x(t), x(t+l), 
x(t+2), x(t+3), x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients Xr(k, t+1) - jXiCk, t+1) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients 
Xr(k, t + 1) - jXi{k, t + 1) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the complex Fourier coefficient 
Xr(k, t) - jXi{k, t) storedin the secondmemory means temporarily , 
obtaining complex Fourier coefficients Xr(k, t+1) - jXi(k, 
t+1) at the frequency interval given by the minimum frequency 
f 1 and the maximum frequency f 2 to the data stream supplied 
since time t+1 based on a transfer function expressed in a 
following equation: 

'[ l-2cos[2jrp]z-' +z-^ J 

where A is a positive constant for providing [x(t+N) - x(t)] 
with an amplitude value, and 

fs[ N-1 J 

11* A recursive discrete Fourier transformation method 
according to claim 10 wherein the positive constant value A 
for providing with an amplitude corresponding to a difference 
between the x(t+N) and the x(t) is capable of being set 
selectively with an inverse number of square root of 1, N or 
1/N. 
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12 . A recursive discrete Fourier transformation method 
according to claim 10 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the 

5 Fourier transformation to the set value k is repeated at every 
time having an equal interval. 

13. A recursive discrete Fourier transformation method 
for obtaining complex Fourier coefficients, wherein data 

10 values Xr{t) + jXi{t), Xr(t+1) + jXi(t + l), Xr(t + 2) + jXi(t+2), 
Xr(t + 3) + jXi(t+3), Xr(t+N-1) + 3Xi(t+N-l), Xr(t+N) + 

jXi{t+N) sampled at times t, t + 1, t+2, t+3, t+N-1, t+N 

(N is a positive integer which is 1 or more) each having an 
equal interval are supplied and with such N data values supplied 

15 since time t as a data stream, a frequency component, which 
is degree k (k is 0 or a positive integer smaller thanN) obtained 
by carrying out complex Fourier transformation on the data 
stream, is obtained such that a real part Xr(k, t) thereof 
and an imaginary part Xi{k, t) thereof are complex Fourier 

20 coefficients, the method comprising: 

a first step of storing the data stream Xx{t) + 3Xi(t) , 

Xj:(t+1) + 3Xi(t + l) , Xr(t + 2) + jXi(t + 2) , Xr(t + 3) + 3Xi(t + 3) , . . . , 

Xr(t+N-1) -I- 3Xi(t+N-l) supplied since time t at time t+N-1 

into a first memory means temporarily; 
25 a second step of obtaining the complex Fourier 

coefficients Xr(k, t) and Xi(k, t) of the data stream stored 

temporarily at the first step; 

a third step of storing the complex Fourier coefficients 

Xr(k, t) and Xi(k, t) obtained at the second step into the 
30 second memory means temporarily; and 

af ourth step of by using the data value Xr( t+N) + JXi(t+N) 

supplied at time t+N, the data value Xr(t) + jXi(t) stored 

in the first memory means temporarily and the complex Fourier 

coefficient Xr(k, t) and Xi(k, t) stored in the second memory 
35 means temporarily and used recursively, obtaining complex 

Fourier coefficients Xr(k, t+1) and Xi(k, t+1) to the data 
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stream supplied since time t+1 with respect to a positive 
constant A for giving an amplitude value corresponding to the 
difference between the Xr(t+N) and the Xr(t) according to a 
following equation: 



5 X,(y!:,r+l) = 



sm 



cos 



sm 



N 



10 14 • A recursive discrete Fourier transformation method 

according to claim 13 wherein by using data value yr( t ) + jYiC t ) 
supplied at time t, data value yr(t+N) + jyi(t+N) supplied 
at time t+N and complex inverse discrete Fourier coefficients 
Yr(k, t) and Yi(k, t) with the real part and the imaginary 

15 part obtained with respect to N data values yr(t) + 3yi(t), 
yr{t+l) + jyi(t+l), yr(t+N-l) + Jyi(t+N-1) supplied from 

time t to time t+N-1, 

a real part Yr(k, t+1) and an imaginary part Yi(k, t+1) 
of each complex inverse discrete Fourier coefficient of N data 

20 values supplied since time t+1 are obtained with respect to 
a positive constant value B for giving an amplitude 
corresponding to a difference between the yr( t+N) and the yr( t ) 
as inverse discrete Fourier transformation data, according 
to following equations. 



25 YXk,t+l)=\^,{kj) + ^\yr{t + N)-yM 



COS 



2^ 
N 
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\Y,{k,t)^^[y{t+N)-y,{t)^sm 



Y, {k,t + 1) = {k,t)+:^M^N)- y, )]| 



cos 



-i^^{k,t)+^M + N)-y,{t)^^sm 



and after the obtained inverse discrete Fourier transformation 
5 data is supplied, discrete Fourier transformation on the 
supplied inverse discrete Fourier transformation data is 
carried out, 

the recursive discrete Fourier transformation 
conducting discrete Fourier transformation by using the 
10 positive constant value A corresponding to the positive 
constant value B. 

15 . A recursive discrete Fourier transformation method 
according to claim 13 wherein k indicating the degree of the 
15 complex Fourier coefficient is set to a desired value and the 
Fourier transformation to the set value k is repeated at every 
time having an equal interval. 



16 . A recursive discrete Fourier transformation method 
20 for obtaining complex Fourier coefficients, wherein with 
sampling frequency as fs, data values Xr(t) + Jxi(t), Xr{t+1) 

+ jXi(t+l), Xr{t + 2) + JXi(t + 2), Xr{t + 3) + 3Xi(t + 3), 

Xr(t+N-1) + jXi(t+N-l), Xr(t+N) + jXi(t+N) Sampled at times 
t, t+1, t+2, t+3, t+N-1, t+N (N is a positive integer 

25 which is 1 or more) each having an equal interval are supplied 
and with such N complex data values supplied since time t as 
a data stream, a frequency interval given with a minimum 
frequency f 1 and a maximum frequency f2 to the data stream 
is regarded as a specified frequency interval, a frequency 

30 interval obtained by dividing the specified frequency interval 
by the N is assumed to be a minimum frequency interval and 
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a result of the Fourier transformation provided by carrying 
out complex Fourier transformation at every minimum frequency 
interval is obtained as a frequency component which is k (k 
is 0 or a positive integer smaller than N) times the minimum 
frequency interval, such that a real part Xr(k, t) thereof 
and an imaginary part Xi(k, t) thereof are complex Fourier 
coefficients, the method comprising: 

a first step of storing the data stream Xr(t) + jXi(t) , 

Xr(t + 1) + jXi{t + l) , Xr(t + 2) + 3Xi(t + 2) , Xr(t + 3) + jXi{t + 3) , .... 

Xr{t+N-1) + 3Xi(t+N-l), Xr(t+N) + jXi(t+N) supplied since time 
t at time t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients Xr(k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients 
Xr(k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the data value Xr( t+N) + 3Xi{ t+N) 
supplied at time t+N, the data value Xr(t) + 3Xi(t) stored 
temporarily in the first memory means and complex Fourier 
coefficients Xr(k, t) and Xi(k, t) stored temporarily in the 
second memory means and used recursively, obtaining complex 
Fourier coefficient XrCk, t+1) andXi(k, t+1) in the specified 
frequency interval given by the minimum frequency f 1 and the 
maximum frequency f 2 to the data stream supplied since time 
t+1 with respect to a positive constant A for giving an amplitude 
value corresponding to a difference between the x(t+N) and 
the x(t) according to a following equation: 



X, {k,t + 1) = |z, {k,t) + ^[x^ {t^N}- x^ (f )]| X cos 



{f2-n)k 



N-1 



+ /1 



X, {k,tU^[x, {t + N)-x, {t)t X sin 



if2-n)k 



N-l 



+ /1 



Z, {k,t + 1) = |z, {k,t)+^[x, {t + n)-x^ (/)]| X cos |2^ 



{f2-fl)k 



N-l 



+ /1 



{f2-fl)k 
iV-1 



+ /1 



17 . A recursive discrete Fourier transformation method 
according to claim 16 wherein by using data value yr(t) + jyiCt) 
supplied at time t, data value yr(t+N) + 3yi(t+N) supplied 
at time t+N and complex inverse discrete Fourier coefficients 
Yr(kr t) and YiCk, t) with the real part and the imaginary 
part obtained with respect to N data values yr(t) + jyi(t), 
yr(t+l) + jyi(t + l), yr(t+N-l) + 3yi(t+N-l) supplied from 

time t to time t+N-1, 

a real part Yr(k, t + 1) and an imaginary part Yi(k, t+1) 
of each complex inverse discrete Fourier coefficient of N data 
values supplied since time t+1 are obtained with respect to 
a positive constant value B for giving an amplitude 
corresponding to a difference between the yr{t+N) and the yr( t ) 
as inverse discrete Fourier transformation data, according 
to following equations. 



COS 



' N 



Y,{k,t) + ^[y{t + N)-y,{t)l^sm 



Y, {k,t +l)^^,{k,t) + ^M+N)- 



COS 



B 



N 



and after the obtained inverse discrete Fourier transformation 
data is supplied, discrete Fourier transformation on the 
supplied inverse discrete Fourier transformation data is 
carried out, 

the recursive discrete Fourier transformation 
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conducting discrete Fourier transformation by using the 
positive constant value A corresponding to the positive 
constant value B. 



10 



18. A recursive discrete Fourier transformation method 
according to claim 16 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the 
Fourier transformation to the set value k is repeated at every 
time having an equal interval. 



19 . A recursive inverse discrete Fourier transformation 
method for obtaining complex inverse Fourier coefficients 
wherein complex data values Xr(t) + jXi(t), Xr(t + 1) + jXi(t+l), 
Xr(t + 2) + jXi(t + 2), Xr(t+3) + 3Xi(t + 3), Xr(t+N-1) + 

15 3Xi(t+N-l), Xr(t+N) + jXi(t+N) sampled at times t, t+1, t+2, 
t + 3, t+N-1, t+N (N is a positive integer which is 1 or 

more) each having an equal interval are supplied and with such 
N data values supplied since time t as a data stream, an inverse 
complex Fourier transformation component, which is degree k 

20 (kis 0 orapositive integer smaller thanN) obtained by carrying 
out inverse complex Fourier transformation on the data stream, 
is obtained such that a real part Xr(k, t) thereof and an 
imaginary part Xi(k, t) thereof are inverse complex Fourier 
coefficients, the method comprising: 

25 a first step of storing the data stream Xr(t)+3Xi( t) , 

Xr( t + 1) + jXi{t + l) , Xr(t + 2) + 3Xi(t-H2) , Xr(t + 3) + 3Xi(t + 3) , - • . , 

Xr(t+N-1) + 3Xi(t+N-l) supplied since time t at time t+N-l 

into a first memory means temporarily; 

a second step of obtaining the inverse complex Fourier 
30 coefficients Xr{k, t) and Xi{k, t) of the data stream stored 

temporarily at the first step; 

a third step of storing the inverse complex Fourier 

coefficients Xr(k, t) and Xi(k, t) obtained at the second step 

into the second memory means temporarily; and 
35 a fourth step of by using the data value Xr ( t+N) + jxi { t+N ) 

supplied at time t+N, the data value Xr(t) + JXi(t) stored 
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in the first memory means temporarily and the inverse complex 
Fourier coefficient Xr(k, t) andXi(k, t) stored in the second 
memory means temporarily, obtaining inverse complex Fourier 
coefficients Xr(]c, t+1) and Xi(k, t+1) to the data stream 
supplied since time t+1 with respect to a positive constant 
B for giving an amplitude value corresponding to the difference 
between the Xr(t+N) and the Xr(t) according to a following 
equation : 



X, [k,t + 1) = |x, {k,t) + ^[xXt+N)-x, 



cos 



10 



X, {k,t) + ^[x, {t + iV)- sin 



X, ik,t + 1) = jz, ik,t) + j[x, {t + N)-x. cos 





2^' 


.sin 









20 . A recursive inverse discrete Fourier transformation 
15 method according to claim 19 wherein k indicating the degree 
of the complex Fourier coefficient is set to a desired value 
and the inverse Fourier transformation to the set value k is 
repeated at every time having an equal interval. 

20 21 . A recursive inverse discrete Fourier transformation 

method for obtaining complex inverse Fourier coefficients 
wherein with sampling frequency as fs, data values Xr(t) + 

JXi(t), Xr(t-Hl) + jXi(t+l), Xr(t + 2) + jXi(t + 2), Xr(t + 3) + 

jXi{t + 3), Xx(t+N-1) + 3Xi(t+N-l), Xr(t+N) + 3Xi{t+N) 

25 sampled at times t, t + 1, t+2, t+3, t+N-1, t+N (N is a 

positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 
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since time t as a data stream, a frequency interval given with 
a minimum frequency f 1 and a maximum frequency f 2 to the data 
stream is regarded as a specified frequency interval, a 
frequency interval obtained by dividing that specified 
frequency interval by the N is assumed to be a minimum frequency 
interval and a result of the inverse Fourier transformation 
provided by carrying out inverse complex Fourier 
transformation at every minimum frequency interval is obtained 
as a frequency component which is k (k is 0 or a positive integer 
smaller than N) times the minimum frequency interval, in the 
form of a real part Xr(k, t) and an imaginary part Xi(k, t), 
the method comprising: 

a first step of storing the data stream Xr(t) + 3Xi{t) , 

Xr( t + 1) + 3Xi{ t + 1) , Xr{t + 2) + jXi(t + 2) , Xr(t + 3) + 3Xi(t + 3) , . • . , 

Xr{t+N-1) + jXi{t+N-l), Xr(t+N) 4- jxi { t+N ) Supplied since time 
t at time t+N-1 into a first memory means temporarily; 

a second step of obtaining the inverse Fourier 
coefficients Xr(k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the inverse Fourier coefficients 
Xr{k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the data value Xr( t+N) + jXi{ t+N) 
supplied at time t+N, the data value Xr(t) + jXi{t) stored 
temporarily in the first memory means and preceding inverse 
Fourier coefficients Xr(k, t) andXi(k, t) stored temporarily 
in the second memory means , obtaining inverse Fourier 
coef f icient Xr(k, t+1) andXi(k, t + 1) in the specif ied frequency 
interval given by the minimum frequency f 1 and the maximum 
frequency f 2 to the data stream supplied since time t+1 with 
respect to a positive constant B for giving an amplitude value 
corresponding to a difference between the x(t+N) and the x(t) 
according to following equations: 

X, {k,t + 1) - jz, (yt,/) + -[x,. [t^N)~ X, X cos {2— 



N-1 



+ /1 
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■^lxXk,t)-\'~[x{t'\-N)-x-{t)].xsm\2 



7t 



(/2-/l)A 



+ /1 



X, {k,t + 1) = jz, {k,t)+ (r + at) - ^, X cos |2^ 



+ /1 



- |z, (/r,? ) + (/ + JV) - X, )]| X sin |2 



(/2-/i)A: 



+ /1 



22 . A recursive inverse discrete Fourier transformation 
method according to claim 21 wherein k indicating the degree 
of the complex Fourier coefficient is set to a desired value 
and the inverse Fourier transformation to the set value k is 
repeated at every time having an equal interval. 



